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We propose a mechanism of energy relaxation for carriers confined in a non-polar quantum dot 
surrounded by an amorphous polar environment. The carrier transitions are due to their inter- 
action with the oscillating electric field induced by the local vibrations in the surrounding amor- 
phous medium. We demonstrate that this mechanism controls energy relaxation for electrons in 
Si nanocrystals embedded in a SiC>2 matrix, where conventional mechanisms of electron-phonon 
interaction are not efficient. 
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I. INTRODUCTION 



Frohlich mechanism of electron-phonon coupling is of 
key importance for carrier energy relaxation in bulk polar 
semiconductors as well as in heterostructures formed by 
different polar materials [TUMI- One of the examples of 
such heterostructures is provided by semiconductor quan- 
tum dots of polar material in a crystalline polar environ- 
ment Recently carrier dynamics in heterostructures 
combining both polar and non-polar constituents, such as 
Si nanocrystals embedded into a silica glass matrix, has 
attracted much attention 3^0- Therefore, it becomes 
interesting to analyze if the polar environment can in- 
duce Frdhlich-like interaction within a non-polar inclu- 
sion. However, the electric field associated with a wave 
of the LO-type propagating in a polar crystal cannot pen- 
etrate through the interface with a non-polar inclusion, 
since the dielectric permittivity, e(u>), vanishes at the fre- 
quency of the LO-phonon mode to = lulo- On the other 
hand, it is well known that amorphous media possess 
high-frequency vibrational eigenmodes of strongly local- 
ized character I^MIol. In this work we will show that 
the local modes in a polar amorphous media can induce 
an alternating electric field even within a non-polar inclu- 
sion and thus provide an efficient channel for relaxation of 
electronic excitations of the inclusion. Such a process can 
be responsible for electron relaxation in Si nanocrystals 
in Si02- We expect this relaxation channel to be pre- 
dominant as conduction-band electron interaction with 
optical phonons via deformation potential is forbidden in 
silicon by the crystal symmetry I 111. 

The rest of the paper is organized as follows. In Scc.HIl 
we describe a simplified model of a non-polar spherical 
quantum dot surrounded by a polar environment and de- 
rive equations for the time of carrier relaxation induced 
by local vibrations in the environment. In Sec. IHII we es- 
timate the parameters characterizing local vibrations in 
a polar glass. In Sec. IIVI we present a calculation using 



the values of parameters describing Si nanocrystals in a 
Si02 matrix. Sec. [V] is reserved for conclusions. Some 
auxiliary derivations are given in the Appendix. 




FIG. I: Schematic illustration of the Si QD surrounded by 
randomly oriented dipoles in Si02 matrix. The radius-vector 
ro is shown for one of the dipoles po- 



ll. INTERACTION OF CONFINED 
ELECTRONS WITH LOCAL VIBRATIONS OF A 
POLAR GLASS 

Let us consider a spherical semiconductor nanocrystal 
embedded into a polar glass matrix. We will characterize 
the jth local vibration mode in the glass by the oscillating 
dipolc moment 
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where ef' and uf\t) are the charge and the displace- 
ment of the i-th ion participating in the j-th local vi- 
bration. The system under study is sketched in Fig. 1, 
where a random orientation of the dipoles is assumed. 

According to electrostatic approximation (see Ap- 
pendix) the potential inside the nanocrystal of the radius 
R and dielectric constant £j n induced by a point dipole 



source with polarization p^8{r — Tq"), positioned in the 
matrix with dielectric constant e ut, can be represented 
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and Yj' M {n) are the vector spherical harmonics intro- 
duced in Ref . [II- 

Neglecting retardation, one can use the Fermi golden 
rule to calculate the rate of intraband transitions of the 
electron confined within the spherical nanocrystal under 
the influence of the quasistatic electric field induced by 
the local vibrations. At low temperatures^ this rate is 
given by 
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where summation runs over all the local modes of the 
matrix, Ei and Ef are, respectively, the energies of the 
initial and the final \ipf) states of the confined elec- 
tron. Neglecting electron tunneling outside the sphere, 
one can write the transition matrix element as 
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Introducing the vibration distribution function as 
P(r , P ,^P 01 Lfc^p M 
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(where Pq = const) summation over vibrational modes 
in Eq. ([3]) can be reduced to 
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The angular integration can be done with the help of 



dn po dn ra [ Po ■ Yf+\n ro )} ■ [p ■ Yf:+} (n ro )] * = 

-^■PqSjj'Smm'- (6) 



Thus, finally we get 
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III. DISCUSSION 

We need to estimate the mean value of the dipole mo- 
ment corresponding to a given local vibrational mode. 
Let us first introduce characteristic values, u and Mo, for 
the mean square displacement and atomic mass of the 
atoms participating in the vibration, respectively (we ne- 
glect the fact that we deal with different sorts of atoms). 
This enables us to write 
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where N is the number of atoms participating in the 

(7) 

vibration, w p h = 



^ ph . At low temperatures 
NM u 2 uj 2 h ~ ftwp h , 



and, therefore, 



NM LOph 
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The dipole moment characterizing this vibrational mode 
is related to u through 

p ~ ewy , (10) 

where the factor 7 depends on the relative orientations 
of displacements for atoms participating in the local 
modefll 

The other glass parameters entering Eq. ([7]) are the 
concentration of dipolcs, Pq, and vibrational density of 
states, p . The concentration of dipoles is of the order of 
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where a is the characteristic atomic scale. The rough 
estimation for the density of states is p ~ 1/ (fi^nh) ■ Nu- 
merical simulations performed for silica glass [§ll9lll0l show 
that the actual phase volume of the high-frequency local- 
ized vibrations is at least by a factor of 5 less than the 
total phase volume. This factor decreases the density of 
states by about an order of magnitude. 



IV. MODEL CALCULATION 

In this Section we will perform a model calculation of 
the relaxation time due to the coupling of electron con- 
fined in a quantum dot with local vibrations of the quan- 
tum dot environment. The quantum dot is treated as a 
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spherical quantum well with the infinite potential bar- 
rier. We consider an electron from a simple band with 
an isotropic effective mass m* and use the effective mass 
approximation. In this case the electron states are char- 
acterized by the radial quantum number n, the orbital 
angular momentum I, its projection m onto an arbitrary 
axis, and a projection of the electron spin. Neglecting 
for simplicity the electron spin, within the effective mass 
approximation one can characterize the electron states 
only by the envelope wave function 

^nim{r) = W —Yi m (6,<p)- f (12) 

where ji (x) is the spherical Bcssel function, the numbers 
ifini arc found from 



jl(<Pnl) = 0, 



and the energy levels 
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are degenerate over m. [l~5| Fig. shows the dependence 
of first six energy levels on the nanocrystal diameter, 2R, 
for 2R > 4 nm, where the density of states mass in Si, 
in* = 0.33 mo was used (too is the free electron mass). 
Comparison of the energy level positions with those pro- 
vided by more accurate calculations [IMT7I suggests that 
the approximation of infinite barriers is reasonable for 
the dot under consideration. 

A transition between two of the levels of size quan- 
tization with the energies (| 13[) is possible when their 
diameter-dependent difference falls within the spectral 
range, 



120 meV < hiu ph < 160 meV 



(14) 



corresponding to high-frequency local vibrations of the 
silica glass IMsflol Vertical arrows in Fig. [5] indicate these 
allowed regions and the corresponding transition time is 
shown near each arrow. The calculation is carried out for 
the set of QD parameters close to that of Si/SiC^ system: 
£i„ = 3, e ut = 12. Our model is too simple to describe 
the spectral dependence of glass parameters. Thus, we 
take them at a fixed value in the middle of the range (fl4|) : 
hojph = 140 meV. Other parameters used are as follows: 
N = 15, 7 = VW, P = (1/N) • 1.5 x 10 22 cm" 3 and 
p(huj) = l/(5fiw p h). The vibrating mass, Mo, was taken 
as M = (l/3)M Si + (2/3)M where M si and M are the 
silicon and oxygen atomic masses, respectively. Since the 
levels (|13[) are degenerate over momentum projection to, 
we have summed the relaxation rate (|7|) over the final 
states and averaged over initial ones. 

Fig. [3] demonstrates that typical value of the relax- 
ation time is of the order of 1 ns. The transition rate is 
proportional to 1/-R, so it changes with the nanocrystal 
diameter slower than the energy difference proportional 
to 1/i? 2 . 
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FIG. 2: Relaxation time between the first six energy levels 
for different diameters of the nanocrystal. The values of n r 
and / are shown for each curve. The transition time for the 
corresponding levels at given diameter is indicated near the 
vertical arrows. 



V. CONCLUSIONS 

We have shown that interaction of electrons confined in 
a non-polar quantum dot with local vibrations of amor- 
phous polar environment provides an efficient energy re- 
laxation channel for the "hot" electrons. We demon- 
strated that this mechanism controls energy relaxation 
in Si nanocrystals in Si02 matrix with diameters in the 
range of 4 - 8 nm, where the energy relaxation pro- 
cess can proceed through emission of a single phonon. 
In this case the relaxation time is in nanosecond range. 
Energy relaxation in quantum dots with smaller diam- 
eters should be controlled by multiphonon processes. 
We expect nonradiativc relaxation there to be slower 
with the rate comparable to that of radiative intraband 
transitionSi 18 i 19 
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Appendix. 

Let us calculate electric filed induced by a point dipole 
Po positioned at at the point Vq outside the sphere of ra- 
dius R. Dielectric constants inside and outside the sphere 
are £j n and £ out , respectively. We are interested in the 
field inside the sphere. Within electrostatic approxima- 
tion the electric field E can be expressed via the scalar 
potential ip as 



E = — grad ip. 
The Poisson equation for the potential (p reads 
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and the boundary conditions at the sphere surface are with 



p\r=R-0 = <P\r=R+0, 
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In the case of homogeneous medium (e; n = e ut) the so- 
lution of (|A2|) is obvious. 
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In what follows it is convenient to use the basis of vector 
spherical harmonics 
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where 



max(r, r ),r < = min(r, r ),n = r/r,n = 



r /r is useful to present (|A4[) as a sum of scalar spheri- 
cal harmonics Yjm- After calculating the derivatives we 
obtain 
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When £j n =/= e out the scalar potential can be written as 
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Substituting (|A8|) and (|A6j) into the boundary conditions 
[3|) we find the coefficients tjM and tjM- 
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which allows to express the scalar potential as |T]) . 
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